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Abstract — Optimal delay-throughput trade-offs for two- 
dimensional i.i.d mobility models have been established in [23], 
where we showed that the optimal trade-offs can be achieved 
using rate-less codes when the required delay guarantees are 
sufficient large. In this paper, we extend the results to other 
mobility models including two-dimensional hybrid random walk 
model, one-dimensional i.i.d. mobility model and one-dimensional 
hybrid random walk model. We consider both fast mobiles and 
slow mobiles, and establish the optimal delay-throughput trade- 
offs under some conditions. Joint coding-scheduling algorithms 
are also proposed to achieve the optimal trade-offs. 



I. Notations 

The following notations are used throughout this paper, 
given non-negative functions f(n) and g(n): 

(1) /(«) = 0{g(n)) means there exist positive constants c 
and m such that /(«) < cg(n) for all n > m. 

(2) /(«) = D.(g(n)) means there exist positive constants c 
and m such that /(«) > cg(n) for all n > m. Namely, 
g(n) = 0(f(n)). 

(3) f{n) = ®{g{n)) means that both f{n) = £2(g(n)) and 
f(n)=0(g(n)) hold. 

(4) f(n) = o{g{n)) means that \\m n ^ f (n) / g{n) = 0. 

(5) /(«) = a>(g(n)) means that \im n ^g(n)/f(n) = 0. 
Namely, g(n) =o (/(«)). 

II. Introduction 

Delay-throughput trade-offs in mobile ad-hoc networks have 
received much attention since the work of Grossglauser and 
Tse [9], where they showed that the throughput of ad-hoc 
networks can be significantly improved by exploring the node 
mobility. Recently the trade-off was investigated under differ- 
ent mobility models, which include the i.i.d. mobility [16], 
[21], [12], [23], one-dimensional mobility [3], [8], random 
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walk [5], [6], [7], [19], hybrid random walk [19] and Brownian 
motion [13]. 

In [23], we demonstrated that the optimal trade-offs for 
two-dimensional i.i.d. mobility models can be achieved us- 
ing rate-less codes when the required delay guarantees are 
sufficiently large. In this paper, we extend the results to the 
two-dimensional hybrid random walk, one-dimensional i.i.d. 
mobility and one-dimensional hybrid random walk models. 
The two-dimensional i.i.d. mobility studied in [23] only mod- 
els the case where the network topology changes dramatically 
at each time slot. However Markovian mobility dynamics may 
be more realistic. Thus the two-dimensional hybrid random 
walk model was introduced by Sharma et al in [19], where 
the unit square is divided into 1 /S 2 small-squares, and mobiles 
move from the current small-square to one of its eight adjacent 
small-squares at the beginning of each time slot (The detailed 
description of the two-dimensional hybrid random walk model 
is presented in Section [TlTb . Since the distance each mobile 
can move is at most 2^/2 /S at each time slot, we can use 
different values of S to model mobiles with different speeds, so 
this two-dimensional hybrid random walk model can be used 
for a wide range of scenarios. Note that the two-dimensional 
hybrid random walk model is the same as the two-dimensional 
i.i.d. mobility model when 5=1. One might wonder why the 
results in [23] are necessary given the results in this paper. 
The reason is that the Markovian mobility dynamics in this 
paper requires a different set of tools than those in [23] and 
as a result, the trade-off in this paper are applicable only when 
S = o(l). Thus, the results in [23] cannot be recovered from 
the results of this paper. We wish to comment that one of 
the main differences between this paper and [23] is that, the 
i.i.d. mobility assumption in [23] allows us to use Chernoff 
bounds to obtain concentration results. However, the random 
walk and other mobility models in this paper require the use 
of martingale inequalities to establish the navel patterns of the 
mobiles. 

In this paper, we will also study one-dimensional mobility 
models. These models are motivated by certain types of delay- 
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tolerant networks [22], in which a satellite sub-network is used 
to connect local wireless networks outside of the Internet. 
Since the satellites move in fixed orbits, they can be modelled 
as one-dimensional mobilities on a two-dimensional plane. 
Motivated by such a delay-tolerant network, we consider one- 
dimensional mobility model where n nodes move horizon- 
tally and the other n node move vertically. Since the node 
mobility is restricted to one dimension, sources have more 
information about the positions of destinations compared with 
the two-dimensional mobility models. We will see that the 
throughput is improved in this case; for example, under the 
one-dimensional i.i.d. mobility model with fast mobiles, the 
trade-off will be shown to be ®(%/D 2 /«), which is better 
than ©(y/D/n), the trade-off under the two-dimensional i.i.d. 
mobility model with fast mobiles. We also propose joint 
coding-scheduling algorithms which achieve the optimal trade- 
offs. 

Three mobility models are included in this paper, and each 
model will be investigated under both the fast-mobility and 
slow-mobility assumptions. The detailed analysis of the two- 
dimensional hybrid random walk model and one-dimensional 
i.i.d. mobility model will be presented. The results of the one- 
dimensional hybrid random walk model can be obtained using 
the techniques used in the other two models, so the analysis 
is omitted in this paper for brevity. Our main results include 
the followings: 

(1) Two-dimensional hybrid random walk model: 

(i) Under the fast mobility assumption, it is 
shown that the maximum throughput per 
S-D pair is 0{y/D/n) when S = o(l) and 
D = «(|logS|/S 2 ), and Joint Coding-Scheduling 
Algorithm I [23] can achieve the maximum 
throughput when S = o(l) and D is both 
C0(max{(log 2 «)|logS|/S 6 ,^logn}) and 
o(n/log 2 n). 

(ii) Under the slow mobility assumption, it is shown 
that the maximum throughput per S-D pair is 
0(f/D/n) when 5 = 0(1) and D = co(\logS\/S 2 ), 
and Joint Coding-Scheduling Algorithm II can 
achieve the maximum throughput when S = 
o(l) and D is both co((log 2 n)|logS|/S 6 ) and 
o(n/log 2 «). 

(2) One-dimensional i.i.d. mobility model: 

(i) Under the fast mobility assumption, it is shown that 
the maximum throughput per S-D is O (y/D 2 fn\ 
given delay constraint D. Then Joint Coding- 
Scheduling Algorithm III is proposed to achieve 
the maximum throughput when D is both co(^/n) 
and o {y/n/ 3 <yTogn) . 

(ii) Under the slow mobility assumption, it is shown 
that the maximum throughput per S-D pair is 
O ( y/D 2 /nj . Joint Coding-Scheduling Algorithm 
IV is proposed to achieve the maximum throughput 
when D is o U/n/ 'log 2 n) . 

(3) One-dimensional hybrid random walk model: 

(i) Under the fast mobility assumption, it is shown 
that the maximum throughput per S-D pair is 



0(J/D 2 /n) whenS = o(l) and D = fi)(l /S 2 ) 7 and 
Joint Coding-Scheduling Algorithm III can achieve 
the maximum throughput when S = o(l) and D 
is both co(max{(log 2 n)|log5|/5 4 ,iynlog«}) and 

(ii) Under the slow mobility assumption, it is shown 
that the maximum throughput per S-D pair is 
0(^/D 2 /n) whenS = o(l) and D = co(l /S 2 ), and 
Joint Coding-Scheduling Algorithm IV can achieve 
the maximum throughput when S = o(l) and D is 
both co((log 2 n)|logS|/S 4 ) and o(V^/log 2 «). 
Note that the optimal delay-throughput trade-off are estab- 
lished under some conditions on D. When these conditions 
are not met, the trade-off is still unknown in general, though 
a trade-off of the two-dimensional hybrid random walk model 
with slow mobiles has been established under an assumption 
regarding packet replication in [19]. We also would like to 
mention that when the step size of the two-dimensional hybrid 
random walk is 1 / y/n, our two-dimensional hybrid random 
walk model is identical to the random walk model studied in 
[6], [7], where the optimal delay-throughput trade-off has been 
obtained. Our results do not apply to this case since the set of 
allowed values for D becomes empty in that case (see (1) (i) 
above). 

The remainder of the paper is organized as follows: In 
Section [Till we introduce the communication and mobility 
model. Then we analyze the two-dimensional hybrid random 
walk models in Section|IV] and one-dimensional i.i.d. mobility 
models in Section [V] The results of one-dimensional hybrid 
random walk model are presented in Section [VI] Finally, the 
conclusions is given in Section [VTTI 

III. Model 

In this section, we first present the models that we use 
for mobility and wireless interference. Then the definitions 
of delay and throughput are provided. 
Mobile Ad-Hoc Network Model: Consider an ad-hoc net- 
work where wireless mobile nodes are positioned in a unit 
square. Assume that the time is slotted, we study following 
three mobility models in this paper. 
(1) Two-Dimensional Random Walk Model: Consider a 
unit square which is further divided into 1/S 2 squares 
of equal size. Each of the smaller square will be called 
an RW-cell (random walk cell), and indexed by (U x ,U y ) 
where U x ,U y 6 { 1 , . . . , 1 /S} . The unit square is assumed 
to be a torus, i.e., the top and bottom edges are assumed 
to touch each other and similarly the left and right edges 
also are assumed to touch other. A node which is in 
one RW-cell at a time slot moves to one of its eight 
adjacent RW-cells or stays in the same RW-cell in the 
next time-slot with each move being equally likely as in 
Figure Q] Two RW-cells are said to be adjacent if they 
share a common point. The node position within the 
RW-cell is randomly uniformly selected. There are n S- 
D pairs in the network. Each node is both a source and 
a destination. Without loss of generality, we assume that 
the destination of node i is node i+ 1, and the destination 
of node n is node 1. 



Fig. 1. Two-Dimensional Random Walk Model 



(2) One-Dimensional I.I.D. Mobility Model: Our one- 
dimensional i.i.d. mobility model is defined as follows: 

(i) There are 2« nodes in the network. Among them, 
n nodes, named H-nodes, move horizontally; and 
the other n nodes, named V-nodes, move vertically. 

(ii) Using {xi,yij to denote the position of node ;. 
If node i is an H-node, yj is fixed and Xj is a 
value randomly uniformly chosen from [0,1]. We 
also assume that H-nodes are evenly distributed 
vertically, so y,- takes values l/«,2/n, . . ., 1. V- 
nodes have similar properties. 

(iii) Assume that source and destination are the same 
type of nodes. Also assume that node ; is an H- 
node if ; is odd, and a V-node if i is even. Further, 
assume that the destination of node ; is node i+2, 
the destination of node 2n — 1 is node 1 , and the 
destination of node 2n is node 2. 

(iv) The orbit distance of two H(V)-nodes is defined 
to be the vertical (horizontal) distance of the two 
nodes. 

(3) One-Dimensional Random Walk Model: Each orbit is 
divided into 1 /S RW-intervals (random walk interval). 
At each time slot, a node moves into one of two adjacent 
RW-intervals or stays at the current RW-interval (see 
Figure |2|. The node position in the RW-interval is 
randomly, uniformly selected. 

Communication Model: We assume the protocol model intro- 
duced in [10] in this paper. Let dist(/,y) denote the Euclidean 
distance between node ; and node j, and r, to denote the 
transmission radius of node ;. A transmission from node i can 
be successfully received at node j if and only if following two 
conditions hold: 

(i) dist(/,;) <n; 

(ii) dist(k,j) > (1 + A)dist(/,y) for each node k ^ i which 
transmits at the same time, where A is a protocol- 
specified guard-zone to prevent interference. 

We further assume that at each time slot, at most W bits can 
be transmitted in a successful transmission. 
Time-Scale of Mobility: Two time-scales of mobility are 
considered in this paper. 
(1) Fast mobility: The mobility of nodes is at the same 



Fig. 2. One-Dimensional Random Walk Model 

time-scale as the data transmission, so W is a constant 
independent of n and only one-hop transmissions are 
feasible in single time slot. 
(2) Slow mobility: The mobility of nodes is much slower 
than the wireless transmission, so W ^> «. Under this 
assumption, the packet size can be scaled as W/H(n) 
for H{n) = 0{n) to guarantee //(n)-hop transmissions 
are feasible in single time slot. 

Delay and Throughput: We consider hard delay constraints 
in this paper. Given a delay constraint D, a packet is said to 
be successfully delivered if the destination obtains the packet 
within D time slots after it is sent out from the source. 

Let Af[r] denote the number of bits successfully delivered to 
the destination of node ; in time interval [0, T] . A throughput of 
A per S-D pair is said to be feasible under the delay constraint 
D and loss probability constraint e > if there exists no such 
that for any n > no, there exists a coding/routing/scheduling 
algorithm with the property that each bit transmitted by a 
source is received at its destination with probability at least 
1 — e, and 

VmPt(^->X, Vi) =1. (1) 

IV. Two-Dimensional Hybrid Random Walk 
Models 

The optimal delay-throughput trade-offs of the two- 
dimensional i.i.d. mobility model with fast mobiles and slow 
mobiles have been established in [23]. In this section, we first 
first extend the results to two-dimensional hybrid random walk 
models. We will obtain the maximum throughput for D = 
(0 (| log S | /S 2 ) , and then show that the maximum throughput 
can be achieved using the algorithms proposed in [23] under 
some additional constraints on D. 

A. Upper Bound 

The upper bound is established under the following assump- 
tions: 

Assumption 1: Packets destined for different nodes cannot be 
encoded together. 
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Assumption 2: A new coded packet is generated right before 
the packet is sent out. The node generating the coded packet 
does not store the packet in its buffer. 

Assumption 3: Once a node receives a packet (coded or 
uncoded), the packet is not discarded by the node till its 
deadline expires. 

Note that Assumption 1 is the only significant restriction 
imposed on coding/routing/scheduling schemes. Next we in- 
troduce following notations which will be used in our proof. 

. A[r]:A[r]=E? =1 A i [r]. 

• b : Index of a bit stored in the network. Bit b could be 
either a bit of a data packet or a bit of a coded packet. 

• d/, : The destination of bit b. 

• Cb : The node storing bit b. 

• tf, : The time slot at which bit b is generated. 

• Li, : The minimum distance between node db and node Cb 
from time slot % to time slot tb+D — l, i.e., 



min dist(<4,Q,)(f). 

t b <t<t b +D-1 



Theorem 1: Consider the two-dimensional hybrid random 
walk model with step-size 5 = o(l) and delay constraint D = 
co ( | log 5 | /S 2 ), and suppose that Assumption 1-3 hold. We have 
following results: 
(1) For fast mobiles, 

48V2WT 



A^TZ 

(2) For slow mobiles, 
An 2 



yfi{VD+\)>E[A[T]]. 



K 2 / 



1)>E[A[T}}. 



(2) 



(3) 



Proof: Let N 1 ™ denote the number of time slots, from 
tb + 1 to tb +D, at which node c/, and db are in the same RW- 
cell or neighboring RW-cells. Then for any L G [0,S/y/n), we 
have 

Pr(L b <L) 
= £ Pr (L h < L\N™ = K) Pr (Af = K) 

K=l 



< 




%L 2 



-K) 



< 



1-1- 



7tL 2 



where the first inequality follows from the fact that the node 
position within a RW-cell is randomly uniformly selected, and 
the last inequality follows from the Jensen's inequality. 

Next we consider E [N l b w ] . Let (Uf(t),Uf(t)) denote the RW- 
cell in which node ; is at time slot t, and (Vf (t) ,Vj' (t)) denote 
the displacement of node ; at time slot t, i.e., 



1, 

0, 



w.p. 
w.p. 
w.p. 



and Vf (f) 



w.p. 
w.p. 
w.p. 



It is easy to see that 

um = 
um = 



f/f(0)+£\f(m) mod- 

m= 1 / 

\ 1 

^(0) + L Vf (m) mod - 



1; 



Further, let (Uf_j(t),Uj_j(t)) denote the relative position of 
node i from node j, i.e., 



\ m= 1 

Ulj(t) = ( 1/^(0) +£t^ y (m) 



mod 



S' 



mod -, 



where 



W-j(t)=vm-v?(t)=< 



and 



VZ-j{t) = V?b)-V]{t) = < 



2. 


w.p. 


1, 


w.p. 


0. 


w.p. 


-1, 


w.p. 


-2, 


w.p. 


2, 


w.p. 


lj 


w.p. 


0, 


w.p. 


-1, 


w.p. 


-2, 


w.p. 



So (Uf_j(t),Uf_j(t)) is the consequence of random walk 
(Vf_j(m),Vl_j(m)) with initial position 

(ujLj(p),uLj(p)) = (um-um> u i(°)- u j(°))- 

Note that node Cb and node db are in the same RW- 
cell if [U x Cb _ db {t),U y Cb _ di (0) = (0,0), and in neighboring RW- 
cells if {U x Cb _ db {t\U y Cb _ db {t)) € {(0,1) J (1,0),(1,1),(0,1/S- 
1), (1/5- 1,0), (1/5- 1,1/5— 1)}. Similar to the argument 
in Lemma [T2l provided in Appendix B, we can conclude that 
forD = co(|log5|/5 2 ), 



qq 

E \Nr\ < ^-5 2 D, 
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which implies that 



Pr(L b <L) < 1-M- — 
36L 2 D. 



99 & n 
2 \ TTF 



(4) 



Based on inequality (@), the proof of inequality (01 is similar 
to the proof of Theorem 3 of [23], and the proof of ^ is 
similar to the proof of Theorem 6 of [23]. ■ 

B. Joint Coding-Scheduling Algorithms 

From Theorem [T| we can see that the optimal delay- 
throughput trade-offs of the two-dimensional hybrid random 
walk models are similar to the ones of the two-dimensional 
i.i.d. mobility models [23]. It motivates us to consider the 
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algorithms proposed in [23]. As in [23], we define and 
categorize packets into four different types. 

• Data packets: There are the uncoded data packets that 
have to be transmitted by the sources and received by the 
destinations. 

• Coded packets: Packets generated by Raptor codes. We 
let (i,k) denote the k th coded packet of node i. 

• Duplicate packets: Each coded packet could be broadcast 
to other nodes to generate multiple copies, called dupli- 
cate packets. We let (i,k,j) denote a copy of (i,k) carried 
by node j, and (i,k,J) to denote the set of all copies of 
coded packet (i,k). 

• Deliverable packets: Duplicate packets that happen to be 
within distance L from their destinations. 

We will show that the optimal trade-offs can be achieved using 
Joint Coding-Scheduling Algorithm I and II presented in [23] 
with the following modifications: 

(1) For the fast mobility model, we use Joint Coding- 
Scheduling Algorithm I with the following modification: 
2D j (25 M) data packets are coded into D/M coded 
packets; 

(2) For the slow mobility model, we use Joint Coding- 
Scheduling Algorithm II with the following modifica- 
tions: D/l data packets are coded to D coded packets. 

For the detail of the algorithms, please refer to [23]. 

Theorem 2: Consider the two-dimensional hybrid random 
walk models. 

(1) Fast mobility model: Suppose that S is o(l), D is both 
Co(max{log 2 rc|logS|/(S 6 ),^log«}) and o(«/(log 2 «)), 
and the delay constraint is 6D. Then under the fast 
mobility model, given any e there exists no such that for 
any n > «o, every data packet sent out can be recovered 
at the destination with probability at least 1 — e, and 



lim Pr 



> 



(—) 

V252(V 




1. 



(5) 



by using the modified Joint Coding-Scheduling Algo- 
rithm I. 

(2) Slow mobility model: Suppose that S is o(l) and D 
is both o(log 2 «|log5|/(5 6 )) and <?(n/(log 3 «)), and the 
delay constraints is 16D. Then under the slow mobility 
model, given any e there exists «o such that for any 
n > no, every data packet sent out can be recovered at 
the destination with probability at least 1 — e, and 



1. (6) 



by using the modified Joint Coding-Scheduling Algo- 
rithm II. 

Proof: Let A denote the area of a cell, and M[t] to denote 
the number of nodes in the cell at time slot t. A cell is said 
to be a good cell at time t if 

9 - - r , 11 - 

-An + l<M[t]<-An. 

Proof of (1): We consider one super time slot which con- 
sists of 6D time slots, and calculate the probability that the 




2D j (25M) data packets from node i are fully recovered at the 
destination, where M = y/n/D is the mean number of nodes 
in each cell. The proof will show the following events happen 
with high probability. 

Node distribution: All cells are good during the entire super- 
time-slot with high probability. Letting denote this event, 
we will show 



Pr( 



(7) 



Broadcasting: At least 16D/(25M) coded packets from a 
source are successfully duplicated after the broadcasting step 
with high probability, where a coded packet is said to be 
successfully duplicated if the packet is in at least 4M/5 distinct 
relay nodes. Letting Aj denote the number of coded packets 
which are successfully duplicated in a super time slot, we will 
first show that 



Pr A: > 



16D 



25M 



> 1 



55D 

n 



p 

- e wm , 



(8) 



Receiving: At least 3D/(25M) distinct coded packets from a 
source are delivered to its destination after the receiving step 
with high probability. Letting B, denote the number of distinct 
coded packets delivered to destination i+ 1 in a super time slot, 
we will show 



Pr Bi > 



3 D 
25M 



16 D 

> > 1 

~ 25 M 



2e 



logD 



' SOOMlogD 



(9) 



From inequalities (|7]), ([8]) and (0, we can conclude that 
under the Joint Coding-Scheduling Algorithm I, at each super 
time slot, the 2D/(25M) data packets can be successfully 
recovered with probability at least 



1 



1 



55£> 



D 

hum/ 



logD 

- 2e 85oo 



D 

500MlogD 



The rest of the proof is the same as the proof of Theorem 4 
in [23]. 

Analysis of node distribution: Since D = o(«/log 2 «) 
implies M = ©(log 77). Inequality (0 can be obtained from 
the Chernoff bound and union bound. 

Analysis of broadcasting step: Consider the broadcasting 
step. Note that when occurs, node i is selected to broadcast 
with probability at least 10/(1 1M) at each time slot. Let 38 ft] 
denote the event that node ; is selected to broadcast in time 
slot t. From the Chernoff bound, we have 



Pr 



> 



9 D 

Tim 



> 1 



D 
600A/ 



(10) 



So node ; broadcasts 9D/(11M) coded packets with a high 
probability. Each coded packet is broadcast to 9M/10 relay 
nodes. 

According to Step (2)(ii) of Joint Coding-Scheduling Algo- 
rithm I [23], each relay node keeps at most one packet for each 
source. Consider duplicate packet (i,k,j). It could be dropped 
if node j is in the same cell as node ; and node i is selected 
to broadcast. Thus, the probability that (i,k,j) is dropped is 
at most 

11 DM 10 1 _ 11 D 
To" X ~9M ~ ~9~^ 



(11) 
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due to the following two facts: 
(a) Let Jffji[D] denote the event that node j is in the same 
cell as node ; in at least one of D consecutive time slots. 
Similar to ©, it can be shown that 



(12) 



under the delay constraint given in the theorem, 
(b) When Sf occurs, node i is selected to broadcast with 

probability at most 10/ (9M) at each time slot. 
Now suppose source i broadcasts D, coded packets, so 
9MD//10 duplicate copies are generated. Let Nf denote the 
number of duplicate packets of node i dropped in the broad- 
casting step. From the Markov inequality and inequality (fTTT i. 
we have 



Pr 



rw MDi 

Nf > 

' ~ 50 



< 



< 



Nf 



50 



(=1 

£f=i h 



Di 



MP, 
50 



9MDj UD 
10 X 9n 



MDi 
50 



55D 



which implies 




1,^=5, >1 



55D 



t=\ 



since otherwise, more than MD, /50 duplicate packets would 
be dropped. Inequality ([8]l follows from inequalities (TOT l and 
(HOll. 

Analysis of receiving step: We group every 3D/logD time 
slots into big time slots, named as b-time-slot and indexed by 
tb, and then divided every b-time-slot into three equal parts, 
indexed by tbi, tbi and ^ 3 as in Figure [3] We first calculate 



One Big Time Slot 



f M '6,2 *b,3- 



One Super Time Slot 



Fig. 3. The Division of A Super Time Slot 



the probability that coded packet (i,k) is delivered in tb2- Let 
i,k) [ f w] denote the event that at least one copy of packet 
(i,k) becomes deliverable in tb.2- If {i>k) is m at l east 4M/5 
relay nodes, we can obtain 



Pr(^ a) [t b , 2 }) 



> 



> 




(14) 



51ogD 
due to the following facts: 
(a) Given D/logD = co(logn| logS\/S 6 ), from Lemma [HI 
provided in Appendix B, we know that with probability 



at least 1 — 1/n, two nodes are in the same RW-cell for 

at least 4DS 2 /(51ogZ>) time slots, 
(b) Given two nodes are in the same RW-cell, the probability 

that they are in the same cell is M/(nS 2 ). 
Next note that the duration of tbi and ^3 are of a larger 
order than the mixing time of the random walk (the mixing 
time is defined in Appendix B). From the definition of the 
mixing time, we have that at any time slot belonging to tb,i, the 
nodes are almost uniformly distributed in the unit square. Let 
&(i.k) [tb,2\ denote the event that coded packet (z, k) is delivered 
to its destination in t/,2- Following the argument used to prove 
inequality (13) of Theorem 4 in [23], we have 



Now let x, denote the positions of the nodes at time slot t, 
and 

X = {Xf} 3(*-l)D . 1 
' iogD +1 

for k — 1, . . . ,51ogD/3. Also let ^uja denote the event 
that (i,k) is delivered in the receiving step. It is easy to 
see that ja occurs if k \ [t^] occurs for some tb .2 € 
{1, . . . ,51ogD/3}. Note that {@(i t k)[tb,2]}t b2 are mutually in- 
dependent given X, so from inequality ( fTBI l, we have 



Pr! 



X) > 1 



Further since B ( - > Y^L\ l@ ( , t) i we can conclude that 



1 



51og£>/3 

2oiog£>y -5' 



16 D 
25M 



> 



16 D 
125 M' 



(16) 



We next bound the number of distinct coded packets deliv- 
erable in tb- Similar to inequality (Tl4"i i. we have 



Pr(^)M) 



< 



IogD 



Note that no two duplicate packets from node ; are in one 
relay node, so [tb] }k are mutually independent. From 

the Chernoff bound, we have 



(DIM 16 
V k— 1 



D 



5 MlogD 



> 1 



D 

400MlogD 



Let &i denote the event that node i obtains no more than 
\6D/(5M\ogD) coded packets at each b-time-slot in the 
receiving step. From the union bound, we have that for 
sufficiently large n, 

'5, 



Pr J 5 / > 1 



■IogD 



>l—e SOOMlogD _ (17) 



Now let Bj(X,Ai,^j) denote the number of distinct coded 
packets delivered to the destination of node i given (X,A,-,^), 
and X tb denote an n x (3D/logD) matrix where the (;,?) entry 
is the position of node i at the f th time slot of b-time-slot tb- 
It is easy to see that the value of Z?,-(X,A,-,<^i) is determined 
by {X rfc }, i.e., there exists a function fi^A-S) sucn mat 

Bj(X,Aj )t ^/) = /(x^JT;) ( X li ■ ■ ■ ! X 51ogD/3) • 



7 



From the definition of 
following condition, 



function /, 



satisfies the (2) For slow mobiles, 



WO 



.X 



fb-i 



^-1 



iYf s ,X fi)+ i, 



^51ogD/3j 



L 51ogD/3; 



4wr A 4 /^ 3 ^( 



D+l)>E[A[T]\. 



< 



16 

5 Mlog 



^18) rrooj: Recall that L b is the minimum distance between 
node di, and node c b from time slot t b to 



It is easy to see that {X. tb } are mutually independent 
given (X,A,-,J^j). Then invoking Azuma-Hoeffding inequality 
provided in Appendix A, we can conclude that 

125 M) 



Pr(2? i (X J A, J ^i)>^[B/|X»Ai,^ 1 ] 
> 1 



logD 

- 2e 



(19) 



holds for any X and A, . Inequality (|9) follows from inequalities 
(H6), (ED and (H9). 

Proof of (2): We consider one super time slot which consists 
of 16D time slots, and calculate the probability that the Djl 
data packets from node i are fully recovered at the destination. 
Let M\ = Un/D, which is the mean number of nodes in each 
cell at the broadcasting step. Following the analysis above, 
we can prove that the following events happen with high 
probability. 

Node distribution: All cells are good during the entire super- 
time-slot with high probability, i.e., 

Pr(Sf) > 1-4- 
n z 

Broadcasting: At least AD/ 5 coded packets from a source are 
successfully duplicated after the broadcasting step with high 
probability. Specifically, we have 



Pr 



4 

A,>-D 



, 50 

> 1 9> 

mV 



(20) 



where a coded packet is said to be successfully duplicated if 
it is in AM\/5 distinct relay nodes. 

Receiving: At least D/6 distinct coded packets from a source 
are delivered to its destination after the receiving step with 
high probability. Specifically, we have 



Pr B 



D 

> — 
~ 6 



4 

5 D 



>l-2e 



(21) 



From inequalities (|7), d20l ) and ( 1211 ), we can conclude 
that under the modified Joint Coding-Scheduling Algorithm 
II, at each super time slot, the D/7 source packets can be 
successfully recovered with probability at least 



— =r — 2e 



logfl 
"TO" 



and theorem holds. 



V. One-Dimensional I.I. D. Mobility Models 
A. Upper Bounds 

Theorem 3: Consider the one-dimensional i.i.d. mobility 
model, and assume that Assumption 1-3 hold. We have fol- 
lowing results: 
(1) For fast mobiles, 



SWT 



TZA 2 



^( 3/ ^D+l)>E[A[T] 



Proof: 

D - 1 . If the orbits 
of node Cf, and dj, are vertical to each other, then Lt,<L holds 
only if at some time slot f, node Cj, and dj, are in the square 
with side length 2L as in Figure [4] In this case, we have 

Vx{L b <L) < 1 - (1 - 4L 2 ) D . 

If the orbits of node c/, and db are parallel to each other, then 
it is easy to verify that 

Pr (L b <L) <1-(1-2L) D . 

Thus for L < 1 /2, we can conclude that 

Pr (L b < L) < 1 - (1 - 2L) D < 2LD. 

The rest of the proof is similar to Theorem 3 and Theorem 6 
in [23]. ■ 



O node i 



node j 



2L 



Fig. 4. The Two Orbits are Vertical to Each Other 

B. Joint Coding-Scheduling Algorithm for Fast Mobility 
Choose 

M= 3 



We divide the unit square into y/n/M horizontal rectangles, 
named as H-rectangles; and y/n/M vertical rectangles, named 
as V-rectangles as in Figure [5] A packet is said to be destined 
to a rectangle if the orbit of its destination is contained in the 
rectangle. 

The algorithms for the one-dimensional i.i.d. mobility model 
has four steps. The first step is the Raptor encoding. The 
second step is the broadcasting step. In this step, the H(V)- 
nodes broadcast coded packets to V(H)-nodes. The third step 
is the transporting step, where the V(H)-nodes transport the 
H(V)-packets to the H(V)-rectangles containing the orbits of 
corresponding destinations, and then broadcast packets to the 
H(V)-nodes whose orbits are contained in the rectangles. After 
the third step, all duplicate packets are carried by the nodes 
that move parallel with the destinations and their orbit distance 



8 




Fig. 5. The Unit Square is Divided into H(V)-Rectangles 



is less than y/M/n. The fourth step is the receiving step, where 
the packets are delivered to the destinations. 

Since duplicate copies are generated in both the broadcast- 
ing step and the transporting step. To distinguish them, we 
name the duplicate packets generated at the broadcasting step 
as B-duplicate packets, and the duplicate packets generated 
at the transporting step as T-duplicate packets. Also we say 
a B-duplicate packet is transportable if it is in the rectangle 
containing the orbit of the destination of the packet. 

Consider a cell with area A and use M H ^ V ^ [t] to denote the 
number of H(V)-nodes in the cell. For the one-dimensional 
mobility model, a cell is said to be a good cell at time slot t 
if 

y^An + 1 < M H(V) [t] < ^An. 

Next we present the Joint Coding-Scheduling Algorithm III, 
which achieves the maximum throughput obtain in Theorem|4] 
Note that in the following algorithm, each time slot is further 
divide into C mini-time slots, and each cell is guaranteed to 
be active in at least one of mini-time slot within each time 
slot. 

Joint Coding-Scheduling Algorithm III: The unit square is 
divided into a regular lattice with n/M cells, and the packet 
size is chosen to be W/(2C). We group every ID time slots 
into a super time slot. At each super time slot, the nodes 
transmit packets as follows. 

(1) Raptor Encoding: Each source takes 2D/(35M) data 
packets, and uses the Raptor codes to generate D/M 
coded packets. 

(2) Broadcasting: This step consists of D time slots. At 
each time slot, the nodes execute the following tasks: 

(i) In each good cell, one H-node and one V-node 
are randomly selected. If the selected H(V)-node 
has never been in the current cell before and not 
already transmitted all of its D/M coded packets, 
then it broadcasts a coded packet that was not 
previous transmitted to 9M/10 V(H)-nodes in the 
cell during the mini-time slot allocated to that cell. 

(ii) All nodes check the duplicate packets they have. If 
more than one B-duplicate packets are destined to 



the same rectangle, randomly keep one and drop 
the others. 

(3) Transporting: This step consists of D time slots. At 
each time slot, the nodes do the following: 

(i) Suppose that node j is a V-node ,and carries B- 
duplicate packet (i,k,j). Node j broadcasts (i,k,j) 
to 9M/10 H-nodes in the same cell if following 
conditions hold: (a) Node j is in a good cell; (b) 
B-duplicate packet (i,k,j) is the only transportable 
H-packet in the cell. 

(ii) Each node checks the T-duplicate packets it car- 
ries. If more than one T-duplicate packet has the 
same destination, randomly keep one and drop the 
others. 

(4) Receiving: This step consists of 5D time slots. At each 
time slot, if there are no more than two deliverable 
packets in the cell, the deliverable packets are delivered 
to the destinations with one-hop transmissions. At the 
end of this step, all undelivered packets are dropped. The 
destinations decode the received coded packets using 
Raptor decoding. 

Theorem 4: Consider Joint Coding-Scheduling Algorithm 
III. Suppose D is (o{y/n) and o {y/n/ 3 ^Togn) , and the delay 
constraint is ID. Then given any e > 0, there exists no 
such that for any n > «q, every data packet sent out can be 
recovered at the destination with probability at least 1—6, and 
furthermore 



1. 



Proof: Consider one super time slot and let £f denote the 
event that all cells are good in the super time slot. The proof 
will show the following events happen with high probability. 
Node distribution: All cells are good during the entire super- 
time-slot with high probability. Specifically, it is easy to verify 
that 




Pr(Sf) > 1 -\. 



(22) 



Broadcasting: At least 2D/ (3M) coded packets from a source 
are successfully duplicated after the broadcasting step with 
high probability, where a coded packet is said to be success- 
fully duplicated if it has at least 4M/5 B-duplicate packets. 
Specifically, we will show 



Pr 



2D 

Ai> 

3M 



>i-12. 

M 



(23) 



Transporting: At least 9D/(10M) coded packets from a 
source are successfully transported after the transporting step 
with high probability, where a coded packet is said to be 
successfully transported if it has at least AM/ 5 T-duplicate 
copies. Letting C, denote the number of successfully trans- 
ported packets from node we will show 



Pr 



9 D 

C, > 

70 M 



D 100 

>l-3e~T500M =-. (24) 

M- 
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Receiving: At least 9D/(140M) distinct coded packets from a 
source are delivered to its destination after the receiving step. 
Specifically, we will show 



Pr 



9 D 

Bi> 

~ 140M 



9 D 

C, > 

70 M 



>l-2e 



D 

800020 



(25) 



If the claims above hold, the probability that the D/IQM 
data packets are fully recovered in one super time slot is at 
least 

1 



1 



D 

- 2e 8000/1/ . 



40 _ d 100 

3e 1500M 

.•i 2 M M 2 
The theorem follows from a similar argument provided in 
Theorem 4 of [23]. 

Analysis of broadcasting step: Assume that occurs, then 
at each time slot, node i is selected with probability 10/(1 IM). 
Note that there are \Jn/M cells on the orbit of node i, and 
node ; is uniformly randomly positioned in one of the cells. 
Thus, the number of coded packets broadcast by node ; is equal 
to the number of non-empty bins of following balls-and-bins 
problem. 

Balls-cmd-Bins Problem: Suppose we have \fnjM bins and 
one trash can. At each time slot, we drop a ball. Each bin 
receives the ball with probability 10/(1 1^/nM), and the trash 
can receives the ball with probability 1 — 10/(1 IM). Repeat 
this D times, i.e., D balls are dropped. 
From Lemma [9] provided in Appendix A, we have 



9D 



,k=l 



11M 



> 1 -, 



We say two nodes are competitive with each other if the 
orbits of their destinations are contained in the same rectangle, 
so each node has \/Mn — 1 competitive nodes. Suppose that 
node i is an H-node and node j is a V-node. Let Nfj(t) denote 
the number of node i's competitive nodes in the V-rectangle 
containing the orbit of node j at time slot t. Since nodes 
are uniformly, randomly positioned on their orbits, from the 
Chernoff bound, we have 



> 1 



_ M 

e ™. 



(26) 



Now consider B-duplicate packet (i,k,j) and assume that node 
Z, a competitive of node i, is in the V-rectangle containing the 
orbit of node ;. Then (i,k,j) might be dropped when it is in 
the same cell as node z, and node z is selected to broadcast. 
The probability of this event is at most 

(27) 

From d26| i, d27b and the union bound, we can conclude that 
the probability that (z, k, j) is dropped at time slot t is at most 




_m 11 M 10 m 11 M 

e 3W-I Mx\ — x = e 300 -| \ —, 

10 V n 9M 9 V « 



(28) 



which implies that the probability of (i,k,j) dropped in the 
broadcasting step is at most 

/ M 11 fM\ D M 11 

1_ l_ e -300 \/— <De~™-\ 



Inequality ( 1231 follows from above inequality and the Markov 
inequality. 

Analysis of transporting step: Consider an H-node ;'. Let 
CWjt) denote the number of B-duplicate packets which are 
contained in the V-rectangle where (i,k) broadcast, and are 
destined to the same H-rectangle as node /. Note the following 
facts: 

(a) Each node has \jMn — 1 competitive nodes. 

(b) Each H-node broadcasts at most D/M coded packets. 
The probability that a coded packet broadcast in a 
specific V-rectangle is at most 



D 



D 



M V n y/nM 
(c) Each broadcast generates 9M/10 duplicate copies. 
Thus, from the Chernoff bound, we have that 

Pr (C (;J .) < DM) > 1 - e _ 36nH , 

which implies that for sufficiently large n, 



Pr {C(i,k) < DM V k) > 1 - -^e"™ 7 > 1 



- g 300M . 



(29) 



Let ^nk) denote the event that a B-duplicate packet is 
broadcast at time slot t in the transporting step. If (i,k) 
is successfully duplicated, i.e., there are at least AM/5 B- 
duplicate copies of (i,k), we have 



~ i ^ , ^ AM M I M 



r l 

L (a)+To— 1 



Further, let UhjA denote the event that at least one copy of 
(i,k) is broadcast in the transporting step. Then for sufficiently 
large «, we can obtain that 



Pr DM] 



> 



1- 1 



AM 




DM+9M/10-1 N 



> -. 

- 4 

Let Cf denote the number of distinct coded packets of node i 
broadcast in the transporting step, i.e., 



C 



D/M 

EM 

k=\ 



(/.A) 



Since different coded packets of node i are broadcast in dif- 
ferent V-rectangles, {^^j are mutually independent. From 
the Chernoff bound, we have 



Pr 



' ~ IM 



2D 

A; -3M' C ^' 



< MD Vk) > 1 — 2e-w<™ . (30) 



In the transporting step, a T-duplicate copy will be dropped 
if the node carrying it obtains another packet destined to the 
same destination. Consider a T-duplicate packet (i,k,l) carried 
by node /. Note following facts: 

(a) Coded packets of node i are broadcast in at most D/M 
V-rectangles. 
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(b) Each rectangle contains at most 9M/10 B-duplicate 
copies from node i. 
Thus, the probability of (i,k,l) dropped at time slot t is at 
most 




The node mobility is independent across time, so the proba- 
bility of (i,k,l) dropped in the transporting step is at most 

MV n \ V n J ~ M 2 

Let Nf denote the number of duplicate packets dropped in 
the transporting step. Note that 9MCf/10 T-duplicate packets 
are generated, and each of them has probability 1 jM 2 to be 
dropped. Using the Markov inequality, we have 



, -rd MC; 

Pr Nf > 



100 



< 



which implies 



Pr 



> 1 



90 

W' 

90 
M 2 



(31) 



since otherwise, more than MCf/ 100 duplicate copies are 
dropped. Inequality (l24t follows form inequality d29t-<f3TT>. 

Analysis of receiving step: The proof is similar to the proof 
of inequality (13) in [23]. ■ 

C. Joint Coding-Scheduling Algorithm for Slow Mobility 

In this subsection, we propose an algorithm which achieves 
the delay-throughput trade-off obtained in Theorem [3] First 
choose 



Mi = 

M 2 = M{ 
and scale the packet size to be 

W 



4D 2 



where c s is a constant independent of n as in [23]. Further, 
we divide the unit square into y«/M2 horizontal rectangles, 
named as H-rectangles; and dnfMi vertical rectangles, named 
V-rectangles. 

Joint Coding-Scheduling Algorithm IV: We group every 
14D time slots into a super time slot. At each super time slot, 
the packets are coded and transmitted as follows: 

(1) Raptor Encoding: Each source takes D/50 data pack- 
ets, and uses the Raptor codes to generate D coded 
packets. 

(2) Broadcasting: The unit square is divided into a regular 
lattice with n/M\ cells. This step consists of D time 
slots. At each time slot, the nodes execute the following 
tasks: 

(i) The nodes in good cells take their turns to broad- 
cast. If node i is a H(V)-node and has never been in 



the current cell before, it randomly selects 9Mi/10 
V(H)-nodes and broadcasts a coded packet to them, 
(ii) Each node checks the B-duplicate packets it car- 
ries. If there are multiple B-duplicate packets des- 
tined to the same rectangle, randomly pick one and 
drop the others. 

(3) Transporting: The unit square is divided into a regular 
lattice with n/M\ cells. This step consists of 2D time 
slots. At each time slot, the nodes do the following: 

(i) Suppose node /' carries duplicate packet (i,k,j), 
which is an H-packet. If node j is in a good cell 
and and (i,k,j) is deliverable, node j broadcasts 
the packet to 9Mi/10 H-nodes in the cell. 

(ii) Each node checks the T-duplicate packets it carries, 
if there is more than one T-duplicate packet des- 
tined to the same destination, randomly pick one 
and drop the others. 

(4) Receiving: The unit square is divided into a regular 
lattice with n/M 2 cells. This step consists of 12D time 
slots. At each time slot, the nodes in good cells do the 
following at the mini-time slot allocated to their cells: 

(i) The nodes which contain deliverable packets ran- 
domly pick one deliverable packet and send a 
request to the corresponding destination. 

(ii) For each destination, it accepts only one request. 

(iii) The nodes whose requests are accepted transmit 
the deliverable packets to their destinations using 
the highway algorithm proposed in [4]. 

At the end of this step, all undelivered duplicate packets 
are dropped. Destinations use Raptor decoding to decode 
the received coded packets. 
Theorem 5: Consider Joint Coding-Scheduling Algorithm 
IV. Suppose D is both o(l) and o (-^/n/log 2 ??) , and the delay 
constraint is 14D. Then given any e > 0, there exists no such 
that for any « > «0i every data packet sent out can be recovered 
at the destination with probability 1—6, and furthermore 




Proof: Following the analysis of Theorem |4j we can 
show the following events happen with high probability. 
Node distribution: All cells are good during the entire super- 
time-slot with high probability, i.e., 



Pr( 



n L 



(32) 



Broadcasting: At least 3D/10 coded packets from a source 
are successfully duplicated after the broadcasting step with 
high probability, where a coded packet is said to be success- 
fully duplicated if it has at least Mi/3 B-duplicate packets. 
Specifically, we have 



Pr 



^15° 



n L 



(33) 



Transporting: At least 3D/40 coded packets from a source 
are successfully transported after the transporting step with 
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high probability, where a coded packet is said to be success- 
fully transported if it has at least AM\/5 T-duplicate copies. 
Specifically, we have 



Pr 



> 1 -« 



d 180 

'ism 



logn 



(34) 



Receiving: At least D/4-0 distinct coded packets from a 
source are delivered to its destination after the receiving step. 
Specifically, we have 



Pr I B t > — 
' ~ 40 



,Q > —d) > 1 -2e"A. 
' ' - 40 



(35) 



Thus, the probability that the D/50 data packets are fully 
recovered in one super time slot is at least 



n L 



d 180 _ D 
~3EWT 2e TWin 



log 71 



and theorem holds. 



like to comment that the key to establishing the optimal delay- 
throughput trade-off is to obtain P, j(D,L), the probability 
that node i hits node j in one of D consecutive time slots 
given a hitting distance L. For example, under the two- 
dimensional hybrid random walk model, the upper bound was 
obtained under the condition D = (o{\\ogS\/S 2 ) since it was 
the condition under which we established an upper bound on 
Pij(D,L) (inequality (HJi). Further, the maximum throughput 
was shown to be achievable under a more restrict condition 
D = o((log 2 n)| log5|/5 6 ) since it was the condition under 
which we established a lower bound on Pjj(D,L) (inequality 
(fT-fl) ). Thus, if we can find techniques to compute Ptj(D,L) 
without using the restricts on D, then the delay-throughput 
trade-offs can be characterized more generally. This is a topic 
for future research. 

Acknowledgment: We thank Sichao Yang for useful dis- 
cussions during the course of this work. 



VI. One-Dimensional Hybrid Random Walk Model, 
Fast Mobiles and Slow Mobiles 

In this section, we present the optimal delay-throughput 
trade-offs of the one-dimensional hybrid random walk model. 
The results can be proved following the analysis of the 
one-dimensional i.i.d. mobility and the analysis of the two- 
dimensional hybrid random walk. The details are omitted here 
for brevity. 

Theorem 6: Consider the one-dimensional hybrid random 
walk model and assume that Assumption 1-3 hold. Then for 
S = o(l) and D = co(l /S 2 ), We have following results: 
(1) For fast mobiles, 



24WT 



1 

7TA2 



v^( J VD+l)>£[A[r]]. 



(36) 



When S = o(l) and D is both 
ft) (max{ (log" n) | log S | /S 4 , \fn log «} ) and 
o{^/n/ 3 tylogn), Joint Coding-Scheduling Algorithm 
III can be used to achieve a throughput same as 
except for a constant factor. 
(2) For slow mobiles, 



l2WT v i 3/ ^"(^+ ! ) ^ e {mt}}- 



(37) 



When 5 = o(l) and D is both ft)((log 2 n)|logS|/S 4 ) and 
o{yfn i 'log 2 n), Joint Coding-Scheduling Algorithm IV 
can be used to achieve a throughput same as ( T3~Tb except 
for a constant factor. 

□ 

VII. Conclusion 

In this paper, we investigated the optimal delay-throughput 
trade-off of a mobile ad-hoc network under the two- 
dimensional hybrid random walk, one-dimensional i.i.d. mo- 
bility model and one-dimensional hybrid random walk model. 
The optimal trade-offs have been established under some 
conditions on delay D. When these conditions are not met, 
the optimal trade-offs are still unknown in general. We would 
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Appendix A: Probability Results 

The following lemmas are some basic probability results. 

Lemma 7: Let X\,...,X„ be independent 0—1 random 
variables such that = fi. Then, the following Chernoff 
bounds hold 

Pr(£z,<(l-5)M^ < e- 52 ^ 2 ; (38) 
Pr^X ; >(l + %j < e- 52 ^ 3 . (39) 

Proof: A detailed proof can be found in [15]. ■ 
Lemma 8: Assume we have m bins. At each time, choose 
h bins and drop one ball in each of them. Repeat this n times. 
Using N\ to denote the number of bins containing at least one 
ball, the following inequality holds for sufficiently large «. 

Pr(_Vi<(l-5)m_5i) < 2e- s2 '" fi ^ 3 . (40) 

where p\ = 1 — e~^. 

Proof: Please refer to [23] for a detailed proof. ■ 
Lemma 9: Suppose « balls are independently dropped into 
m bins and one trash can. After a ball is dropped, the 
probability in the trash can is 1 — p, and the probability in 
a specific bin is p/m. Using N2 to denote the number of bins 
containing at least 1 ball, the following inequality holds for 
sufficiently large n. 

Pr(N z <(l-5)mp 2 ) < 2 e - s2m ^ 3 ; (41) 
_ n p 

where p2 = 1 — e «< . 

Proof: Please refer to [23] for a detailed proof. ■ 

Next we introduce the Azuma-Hoeffding inequality. 

Lemma 10: Suppose that Xq, . . . ,X n are independent ran- 
dom variables, and there exists a constant c > such that 
/(X) = f(X\ ,X n ) satisfies the following condition for any 
i and any set of values jci,... ,x n and y, : 

(/"(xi , . . . ,Xj— l ,Xi : Xi-{-l , . . . 

f(x u . . . y,-,x /+ i , . . . ,x„)\ < c. 
Then we have 

Pr(\f(X)-E\f(X)]\>S)<2e 
Proof: A detailed proof can be found in [15]. ■ 



Appendix B: Properties of Random Walk 

Consider following two random walks. 

(1) One dimensional random walk: A random walk on a 
circle with unit length and 1/5 points. At each time 
slot, a node moves to one point left, one point right or 
doesn't move with equal probability as in Figure [6] 




1/3 



Fig. 6. One Dimensional Random Walk 

(2) Two dimensional random walk: A random walk on a 
unit torus with 1/S 2 points. At each time slot, a node 
moves to one of eight neighbors or doesn't move with 
equal probability as in Figure [7] 



• • • • f 




• • • • • 



Fig. 7. Two Dimensional Random Walk 

We introduce following definitions. 

• Transition matrix P : P = [flj] where Pi j is the probabil- 
ity of moving from point i to point j . 

• Stationary distribution II : A vector which satisfies the 
equation IIP = II. 

• Hitting time 7/,(i,j) : Time taken for a node to move from 
point i to point j. 

• Mixing time T m : 

r m =infsu P £|^-nj| <5 4 , 
' i j 

where PL is the (i,j) th entry of P'. 
Lemma 11: For the one dimensional random walk, we have 
. E[7ft(i,j)] = 0(l/5 2 ). 
. 7;, = 0(|log5|/5 2 ). 

For the two dimensional random walk, we have 
. J B[r / ,(i,j)]=0(|log5|/5 2 ). 
. 7;„ = 0(|log5|/5 2 ). 
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Proof: Please refer to [11] for the one dimensional 
random walk. The proof of the hitting time of the two 
dimensional random walk is presented Lemma 13 in [6], and 
the mixing time result holds since the two dimensional random 
walk can be regarded as two independent one dimensional 
random walks. ■ 
Lemma 12: Let N^y^D] denote the number of visits to 
point j in D time slots starting from point i and ending at 
point k. If D = 0)(|logS|/S), we have 



9 ~ 11 ~ 
— DS <E\N iik \D}}< —DS 

10 ~ 1 J ' 1 JJ ~ 10 



(42) 



where S = S for the one dimensional random walk and S = S 2 
for the two dimensional random walk. Furthermore, if D = 
K<x\\ogS\/S 3 where K = ft)(l), then we have 



Pr 



6 „ 4 
-DS > N iiiM [D]\ > -DS ) > 1 - 2e~s55 . 



(43) 



holds for any X. Then from the Azuma-Hoeffding inequality 
(Lemma ITOb, we have that 

9 



Pr 



< — DS 
100 



> 1 -2e~5i5. 



for any X, and inequality d43l holds. 



Proof: First we have 

iVij,k[0]-l 

T h (i,j)+ £ TlQ,])+T h Q,k)=D, 
i=i 

where 7)'(j,j) is the time duration between I th visits to point j 
and (/ + l) th visits to point j. Taking the expectation on both 
sides, we have 

£[r A (iJ)]+£[Mj,k[0]]£[r ft (jJ)]-£[r A (j,j)]+£[7) ! (j J k)]= J D, 
which implies 

prA7 rnll D-E[T h (i,j)]-E[T h (j,k)}+E[T h (U)} 

m - lk[D]] = smm ' 

Inequality (l42l follows from the facts that ^[r/^ij)] = 
0(\logS\/S) wdE[T h (M)] = l/$. 

Next let x[f] denote the position of the node at 

time slot t, X denote {x[f]} for t = 1,^ + 1,2^ + 
l,...,D- 
l,^,...,min 



D, and X m denote {x[t}} for t = 

(m+l)DS : 
a 



mDSf 
a 

Q 



where m = 0,...,a/5~ — 1. 
Further, let x [D] denote the number of visits to point j given 
X. It is easy to see that for any X there exists a function / x 
such that 

^,xM=/x(Xi,...,X a/ ^_ 1 

where {X„,} are mutually independent given X. Note that X,„ 
contains the position information from time slot mDS 2 /a + 1 
to time slot (m+l)DS 2 /a, so 



fx ( x o,-- 
fx ( x o>- 



i X m— 1 1 X m i Xn+ 1 > ■ ■ ■ > X 
• 7 X?7— 1 ; Y m , X m ^ 1 , . . . , X 



a/S 2 -l 



a/S 2 -l 



< 



DS- 
a 



Next note that DS 2 /a = o(|log5|/5), so from inequality 
d42l , we can conclude that for any i, j and k, 



9 DS 3 
which implies that 



<E 





DS 2 










a 





< 



1 1 DS i 

To - ^ - 



—DS<E 

io - 



<^DS 
~ 10 



